On Mittag-Leffler function and associated polynomials 

D. Babusci/'Q G. Dattoli, 2 fl and K. G6rska 3 ' 4 fl 
1 INFN - Laboratori Nazionali di Frascati, 
via E. Fermi, 40, IT 00044 Frascati (Roma), Italy 
2 ENEA - Centro Ricerche Frascati, via E. Fermi, 
45, IT 00044 Frascati (Roma), Italy 
3 H. Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, 
ul.Eljasza-Radzikowskiego 152, PL 31342 Krakow, Poland 
^ENEA - Centro Ricerche Frascati, v. le E. Fermi, 
45, IT 00044 Frascati (Roma), Italy 

Abstract 

The Mittag-Leffler function plays a role of central importance in the theory of fractional deriva- 
tives. In this brief note we discuss the properties of this function and its connection with the 
Wright-Bessel functions and with a new family of associated heat polynomials. 
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The Mittag-Leffler function (MLF) l|, |2| plays a centra 



role in the theory of fractional 
1 

3[ but the increasing interest for 



derivatives (FD) 3 . It has been thoroughly investigated 
FD in applications 4j demands for further studies, eventually leading to further properties 
or to more efficient methods of computation and analysis. 

The importance of FD differential equations to model anomalous diffusion 4|, |5[ has 
been recognized in the past two decades when it has been shown that the evolution of the 

nn 

distributions associated with this process can be modeled by equations of the type p, [7] 

t- c 



d?F a (x,t)=L FP F a (x,t) + 



r i 



a) 



/(*) 



F a (x,0) = f(x), 



(1) 



with a G (0, 1), and Lpp being the Fokker-Plank (FP) operator associated with the specific 



problem under study. The use of the evolution operator method 
an efficient tool for searching the solution of this kind of problems 



ras been proved to be 
8[. Therefore, limiting 



ourselves to the case L FP = k a d%, we can cast the solution of Eq. (JTj) in the form 



F a (x,t) = E a (k a t a d 2 x )f(x 



E a (y) = J2f7Y 



n=0 



a n) 



(2) 



where E a (y) is the MLF. The operator E a (k a t a d%.) is the evolution operator for the FP 
problem in Eq. (pQ) and the associated rules can be found in Refs. (si, 9]. 

Within the framework of FD, the MLF plays the same role of the exponential function 
in the case of ordinary calculus. In the case f(x) = x n , from Eq. (TSJ) we find 



F a \X,tj a FI n (x^k a t j 



[n/2] 
m=0 



X 



n-2m 



(k a t c 



n — 2m)\ T(l + am) ' 



(3) 



where a H n (x, y) is the fractional heat polynomials of order a. Their names derive from the 
so-called heat polynomials (sometimes Hermite heat polynomials), which are the natural 
solutions of the ordinary heat equation. The generating function of this family of polynomials 
is 



y2^- a H n (x,y) = E a (ey) e* 



(4) 



n=0 



and, therefore, they are members of the Appell family 10| satisfying the recurrences 



d xa H n (x,y) = n a if n _i(x, y), 

dy a H n (x, y) = n (n - 1) a H n _ 2 (x, y) 



r(l-a) 



x n . 



(5) 
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To explore more in depth the nature of these polynomials, we consider a function F(x) that 
can be expanded in terms of them (we assume that y plays the role of a parameter) 

oo 

F ( x ) = 5^ a ' n ^ H n{x, y), (6) 

n=0 

with the coefficients a n to be determined. In the case of Eq. flS]), it can be shown that 

f°° 1 t ( t \ 

a H n (x,t a ) = J dsn a (s,t) 1 H n (x,s) n a (s,t) = - gl+1/a g a ( ^> ) > ( 7 ) 



with g a (x) being one-sided Levy stable distributions (see Eqs.(2), (3), and (4) in Ref. 
Therefore, we can write 

POO 

F(x) = dsn a (s,y 1/a )G(x,s) (8) 
Jo 

where 

00 

G(x,s) = iH n (x,s). (9) 

n=0 

We can accordingly conclude that a function can be expanded in series of fractional heat 
polynomials whenever a corresponding expansion in series of ordinary heat polynomials 
exists. The polynomials a H n (x,y) cannot be considered orthogonal, however the identity 
(jHJ) establishes a close links with the ordinary case and with the relevant orthogonality 
properties. Further comments on these aspects will be presented elsewhere. 



jrop 



The function E a (—x 2 ) can be considered as the generalization of Gaussian and its 
erties can be studied using an operational technique recently suggested in Refs. 13 
According to this technique, by defining an operator c such that 

c a (p(0) = ip(a+ 1), (a E R) (10) 

for <f (fx) = l/r(/i) we can write 

00 

E a {~X 2 ) = ^(-l) fc X 2fc C«V(0) = 1 + dax2 m- (11) 
k=0 

Albeit apparently trivial, if correctly used, last equation can usefully be exploited to draw 
many consequences. Let us, for example, consider the integral 



00 



/ dxE a (-x 2 ). (12) 
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As a consequence of Eq. (ITT]) , taking into account the identity 

1 1 



dse- sA s v '\ 



(13) 



and treating c a as a constant, it's easy to show that 

j(2) = . 



7T 



rh-- 



(14) 



The same procedure can be used to evaluate successive derivatives of the function E a (—x 2 
From Eq. fflTl) . one has 

(— l) fc (2 yfc)! 



d n x E a { x 2 ) Y^ { 2 k _ n)[v{l + ak) 
According to our formalism, the function defined by 

^,/3(*)=C Q -V%(0), 



Jlk-n 



X 



(15) 



(16) 



has the following series expansion 

oo 



X 



k=0 



(17) 



i.e., it is the Wright function 



15 



161 ] . By using the representation (TIB"]) for this function, it 



is straightforward to derive the following identities 



d?W a , p {x) = W, 



and 



dxW n r(-x 



7T 



r [a - 



PY 



(18) 



(19) 



Let us now consider the family of polynomials of Laguerre. It is possible to define a 
fractional version of the two-variable type as follows 

^ -Qfc ^.oik y-n—k 



a L n (x,y) = n\ ^ 



k=0 



k\ (n-k)\T(l + ak) 



(20) 



These polynomials satisfy the differential equation 



d y F a (x,y) = a D L>x F a (x,y), 



F a (x,0) 



(-l) n x na 
T(l + na) 



(21) 



where the fractional Laguerre derivative a D Lx = — <9" x d x has been introduced (for 
the notion of the ordinary Laguerre derivative see Refs. (13, lia|). This implies that the 
fractional Laguerre polynomials can be defined by means of the identity 

° L ^= e "" DL "mr™)' (22) 

from which we find 

oo 

a L n (x,y) = e""- 6 *- E a (-x a t). (23) 

n=0 

This generating function can also be obtained by replacing the exponential function with 



15] , and, therefore, 



MLF in the generating function of the ordinary two- variable polynomials 
we can write 

g t%i „ (x ,, )= _*H__) < 24 > 

from which we conclude that the function on the l.h.s. of this equation is the solution of 
Eq. (1211) with initial condition E a (—x a t). Furthermore, it is also easily proved that 

~ { aL n (x, y) = e yt W ha {-x a t). (25) 

n=0 
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